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ABSTRACT 

We prove that static fluid stars can stably support magnetic fields (within the ideal 
MHD approximation). 



1. Introduction 

Convectively unstable stars generate magnetic fields by a dynamo mechanism. When the fluid 
motion stops, should the field die out? 

A necessary condition for the existence of the magnetic field in a static star is the ideal mag- 
netohydrodynamic (MHD) stability of the field. One might think that ideal MHD stability of a 
given equilibrium magnetic field in a static star is an easy problem. The equation for small per- 
turbations near the equilibrium is self-adjoint, and a variati onal principle that decides the stability 
question is easily formulated. But in practice ([Tavleii (|l98d ) and references therein) the variational 
principle has only been used to prove that certain magnetic configuratio ns are unstable. Numerical 
simu lations do suggest that mixed toroidal-poloidal fields can be stable (iBraithwaite and Nordlund 



Here we prove stability of a certain configuration - mixed toroidal-poloidal field in a stably 
stratified incompressible ideally conducting star. The particular configuration that we prove stable 
is artificial - our choice of the configuration was obviously driven by simplicity. However, "by 
continuity", a large class of stable fields close to our con figuration should exist. In fact, the 
numerically stable fields of (jBraithwaite and Nordlundll2005l ) are close to our configuration. 



2. Stable Magnetized Star 
2.1. The initial field 



We will show that a stable equilibrium configuration exists which is close (in the sense specified 
below) to the non-equilibrium initial configuration with the density field p and magnetic field B 
given by 

Pi, r < R 

p{r) = { p2, R<r <Rs (1) 
0, Rs<r 
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B(r) 



^o(^^cos0, -{j[{kr) + =^)sme, ji{kr)sme), r<R 



0, 



r > R 



(2) 



We use spherical coordinatization {r,9,(j)); ji and j[ are the spherical Bessel function and its 
derivative; ji{kR) = 0. 

The magnetic field ([2]) is the field of the minimal energy for a given helicity with a zero normal 
component on the surface of the sphere of radius R, 



V xB = kB. 



(3) 



The importanc e of helicity for the probl em of stability is well known and explained in the next 
subsection (see [Arnold and KhesinI (jl998l ) for discussion and references) . 



The free parameters of the initial configuration are the radius of the star Rg , the radius of the 
magnetized region R, the densities of the magnetized and unmagnetized regions pi and p2, and the 
amplitude of the field Bq. To simplify the proof of stability we assume that the stratification is 
strongly stable, pi 3> P2, and the magnetic field is weak, Gp\R^ ^ Bq. 



2.2. The stable field 

Assume incompressibility and ideal conductivity. Then an equilibrium configuration is an 
extremum of the total (gravitational plus magnetic) energy under incompressible deformations of 
the star and the magnetic field (comoving density constant, magnetic field frozen into the fiuid). 
Stable equilibrium is the minimum of the energy under incompressible deformations. 

The initial configuration is not an equilibrium, because the gravitational energy is at a min- 
imum, while the magnetic energy is not extremal. The Ampere force j x B does vanish inside 
the magnetized region due to ([3]). But there is a singular Ampere force on the boundary of the 
magnetized region at r = i?. This force acts in the direction of deforming the sphere into an oblate 
spheroid. 

Consider the set of all configurations which can be obtained from ()l|2p by incompressible 
deformations. We will show that the minimum-energy configuration belonging to the set is close 
to the initial configuration. In particular, the boundary of the magnetized region of the minimum- 
energy configuration is an axisymmetric oblate spheroid close to a sphere (in the mean square 
deviation sense). 

The idea of the proof is straightforward. A deformation of the star causes a deformation of 
the boundary of the magnetized high-density region. Let ip{9, <j)) be the deviation of the deformed 
boundary from the sphere: 

r = {l + ij{e,(P))R. (4) 
For pi ^ p2, the change of the gravitational energy dWg is a functional of ^ only, after normalizing 
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to the gravitational energy of the initial density ([T]): 

6Wg = \Wg\Wg[llJ\, (S) 

The change of the magnetic energy 5Wm depends on the shape of the boundary ijj and on the 
deformation field inside the magnetized region. But since the frozen-in deformations conserve 
the magnetic helicity, a lower bound on the change of the magnetic energy can be obtained by 
minimizing the magnetic energy for a given boundary and helicity. The boundary condition for 
this minimization is still the vanishing of the normal component of the magnetic field, because this 
property is preserved by the frozen-in deformation of the field. Thus we have 

5Wm > WmWm[lp], (6) 

where Wm is the magnetic energy of the initial field ([2]), and Wmlip] is the dimensionless change of 
the minimal magnetic energy for a given helicity due to the variation of the boundary ip. 

The dimensionless change of the gravitational energy Wg[ip] is positive, and as we show in 
Appendix I, up to the second order in ip, 

Wg[ij] > Cil, 1 = j dnij"^. (7) 

Here and below C denotes positive dimensionless constants. It is assumed that the deformation 
tp does not contain a dipole mode (translation), for this mode clearly does not contribute to the 
magnetic energy, and will be stabilized by the non-zero outer density p2- 

The dimensionless change of the magnetic energy ifmiV'] can be negative, but as we show in 
Appendix II, also up to the second order in ip, 

Wmm>-{C2l'^^ + C3l) (8) 



The total energy of the minimum-energy configuration is smaller than the initial energy, 

dWg + 6Wm < 0. (9) 
It follows that the minimum-energy configuration is close to the sphere: 

J dni^^ < c-^ « 1. (10) 

As shown in Appendix II, the minimum-energy configuration is axisymmetric. The magnetic field 
of the minimum-energy configuration is close to the initial field ([2]). 
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3. Conclusion 



Thus static stars can stably support magnetic fields. Our field is hidden inside the star, but it 
seems likely that a poloidal component which sticks through the surface of the star can be added 
without causing an instability, at least a component with a small enough amplitude. This is because 

the minimum-energy state that we have presented forbids any deformations of the magnetized region 
other than motions along the magnetic surfaces. Therefore, adding a weak net poloidal flux which 
goes through the surface of the star and through the center of the star is likely to cause just a small 
deformation of our equilibrium - one cannot kill this poloidal component by the motions along the 
magnetic surfaces. 



I thank Peter Goldreich, Jeremy Goodman, Andrew MacFadyen and Boris Khesin for discus- 
sions. This work was supported by the David and Lucile Packard foundation. 



A. Gravitational Energy 

Take R = 1, pi = 6, p2 = 0. Then the gravitational potential is 

-3 + r2 + ^ai„yyh„, r<l+V 
\ -2r-i + E&/mr-('+i)l1„, r>l + V, ^ ' 

where a, b are constants, Y are the spherical functions. From continuity of the potential and its 
derivative across the boundary r = 1 -|- '0, we get, to first order in ip, 

-6 

aim = km = 2^ ^ 1 ^'^' (^^) 

where 

= X] ^IrnYlm- (A3) 

Note that ipo is a second order quantity - from volume conservation, to second order in ■0, 

i>i 

To second order in ip, the continuity of the potential gives 

aoYo + ^{aim - bim)Yim = 3V'^. (A5) 
l>i 

which gives 
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The gravitational energy is 

Wg = ^ j d^r p(l) = 3 j^^^ dr sin 6 dO d(l){-3 + + ^ a^^r^;^). (A7) 

The change of energy is, up to second order in V', 

Wg[i;](x j sine de d(t){ilj + 2^)"^ + ^ai^Yi^{j^+ij)^ = ^(l^^^p+^^^a*^)^^^^^ (A8) 

and using ()A21 IA6P , we get the final answer 

^,[V']oc^(l--^)|Vi„^p. (A9) 
i>i 

In the P2 = approximation, there is no penalty for the simple displacement of the high-density 
region (the dipole, / = 1, modes). A non-zero outer density adds a positive contribution from the 
dipole. However, we don't need to include this contribution, because the dipole does not contribute 
to the magnetic energy WmW\- 



B. Magnetic Energy 

Consider magnetic fields that minimize magnetic energy, W = J d^rB'^, for a given helicity, 
H = J d^rA ■ B, where V x A = B. The minimizing fields are solutions of V x B = fcB. The 
minimizing field of the sphere is ([2]). When one deforms the sphere without changing the enclosed 
volume, the value of the Lagrange multiplier k changes, and from W = kH, one gets the new 
minimal energy for a given helicity. 

In our calculation of the minimal energy we prefer to keep k constant by allowing volume 
changing deformations of the boundary. Then, by a uniform expansion or contraction, we restore 
the volume, and calculate the minimum energy from the scaling 

W = CHV-^I'\ (Bl) 

where C depends on the shape of the boundary, but not on the enclosed volume V . 
Thus, we will calculate magnetic field B satisfying 

VxB = B, r <xi + ip{e,(t)), (B2) 

and tangent to the boundary r = xi + ip{0, (j)). Here xi is the first zero of ji, and the unperturbed, 
zeroth order in ■0, solution is 

Bo = f ^ cos 9, - {j[ + ^-^)sme, ji sin o] , (B3) 
V r r J 
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The exact solution for the arbitrary change of the boundary ip (with an appropriate uniform ex- 
pansion harmonic tpo to keep k = 1) is given by ()B2p : 

B = Bo + V Birr, ( l{l + l)-Ylm, ^JlYl^ + {j[ + -)deYirr,, " jldeYirr, + + -)Yir^ . 

^-^ \ r smt) r smfe' r J 

(B4) 

The field must be tangent to the boundary: 

Br = -Bedei} + -^B^d^il:, r = xi + </>), (B5) 
r r sm a 

or, up to second order in ip, 

Z BimKl + l)jiYim = -m - ;^)ii^5e(sin2 9^) + j[^d^i; 

-i, Y: BlrJil + l)j'iYira + 89^ E + 0/ + t)^^'^'™) (B6) 



The argument of Bessel functions here and below is xi. To first order ()B6p gives 

BirrMl + ^)jiYim = "Ji -^^5^ (sin^ 0^), (B7) 



which allows to express the B harmonics in terms of ip harmonics to first order. Using the formula 
for the derivative of spherical functions 

('7_Ln2_ 2\l/2 //2_»„2\l/2 

Sin OdeYirr, = I ( \'^^ _ ^ J - + 1) [-^^ ] Yi-i,m, (B8) 



and integrating ()B7p multiplied by Yj^ over angles we get 

Blm = —Xlm, Xlm ^ T—T ah , -1 ^9 T Vi+l.m " T T V'Z-l.m- (B9) 



The first-order result (|B7p allows to re- write (jB6p . to second order, in the following useful form 



where 



V Bimlil + l)jiYim = -^deFi + -^d^F2, (BIO) 
^ smO sm^^ 



Fi = sin^ e + smedeG2 + F2 = sin^ 0^ - smOdeGi + ^./.Gs)^ (Bll) 

Gl = V BirrrjlY^m, ^2 = V S,^(i; + — (B12) 

Choosing different values for Bi^^ one gets functions with arbitrary / > 1 harmonics but with 
vanishing / = and / = 1 harmonics in the left-hand side of (jBlOp . The right-hand side of (|B10p 
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automatically gives zero for the I = harmonic. The vanishing of the I = 1, m = harmonic means 
that the / = harmonic of Fi vanishes. This gives, to second order in ip, 

V^O = 4f^2,0 + f (— + + l)IX/m|^ + imijI^Xlra) ■ (B13) 

V5 4^7r ^ \ xi Jl J 

We do not need the / = 1, m = ±1 equation, for it simply constrains the orientation of the 
unperturbed field. 

Equation (jBlSp describes the deformations of the sphere, which leave the minimal energy for 
a given helicity unchanged: W = H, because V x B = B. However the deformation (jBlSh changes 
the volume inside the boundary: the volume inside r = xi + 1^(6, (j)) is greater than the volume 
inside the sphere r = xi by 

5V oc 2y^xiV'o + 1^'™!'- (^1^) 

Restoring the volume by a uniform contraction, we get the change of the minimal energy per 
unit helicity caused by the ^-deformation. From (jBl|l : Wm[i^] oc 5V . Using (IB13[) in (IB14jl . we 
obtain the final expression 



JV] OC + Yl (l^'^l^ + ^^l^"^V'LX/m + ^(1 + + • (B15) 



The linear term of this functional, oc V'2,0) says that magnetic energy can be reduced by a 
deformation with a negative (2,0) harmonic, that is by an oblate deformation oc 1 — Scos^^. 

For / ^ 1, the quadratic term in (|B15p is non- negative (minimize over x as if x were indepen- 
dent, using ji{x) oc x^). Therefore, one has a lower bound: it'mlV'] > —(6*2/^''^ + C3/), I = j dO, ifP' . 

Modes with different m are not coupled to each other in (jBlSP : it follows that the minimal 
energy state is axisymmetric. 
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